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RATE-ENHANCED  BURNING  OF  A  PROPELLANT 
OBEYING  A  PYROLYSIS  LAW 

by  F.  A.  Williams 


I.  INTRODUCTION 

It  has  been  found  empirically^  that  by  inserting  wires  in  a 
solid  propellant  grain,  the  burning  rate  of  the  propellant  can  be  in¬ 
creased  by  a  factor  of  5  or  10.  The  shape  of  the  burning  propellant 
surface  becomes  a  series  of  cones,  with  their  vertices  centered  on  the 
wires.  Heat  conduction  from  the  hot  gas,  along  the  highly  thermally 
conductive  wire,  into  the  propellant,  serves  to  maintain  the  conical 
shape  of  the  surface.  The  phenomenon  provides  one  useful  method  for 
obtaining  increased  propellant  burning  rates  in  a  motor  that  is  sub¬ 
jected  to  various  geometrical  restrictions. 

The  purpose  of  the  present  paper  is  to  develop  a  theoretical 
understanding  of  the  mechanism  of  this  method  of  burning  rate  enhance¬ 
ment  and  to  derive  approximate  formulae  that  may  be  used  to  obtain 
estimates  of  the  degree  of  enhancement.  A  theoreticeil  analysis  con¬ 
taining  all  essential  physical  features  of  the  phenomenon  has  been 
given  previously  (Ref.  l)j  but  the  author  believes  that  this  analysis  ' 
is  somewhat  too  conpllcated  and  hot  entirely  perspicuous.  The  present 
study  is  intended  to  clarify  the  nature  of  the  phenomena  analyzed  in 
Ref.  1  and  to  develop  slnpler  working  formulae  that  are.  Inherently, 
Just  as  accurate  as  those  of  Ref.  1. 

In  Section  II  we  shall  consider  a  sinple  problem  related  to 
rate-enhanced  burning  in  order  to  introduce  the  underlying  character 
of  the  process.  The  successful  model  of  rate-enhanced  burning  is  de¬ 
fined  and  discussed  in  Section  III.  In  Section  IV,  what  e^pears  to  be 
the  sinplest  conceivable  ansJ-ysis  of  rate-enhanced  burning  is  presented; 
a  more  conplete  phenomenological  theory,  which  Includes  all  of  the 
physical  effects  included  in  the  analysis  of  Ref.  1  but  which  is  less 
coii;>licated  than  the  analysis  of  Ref.  1,  is  presented  in  Section  V. 
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In  Section  VI,  the  predictions  of  the  theory  of  Section  V  are  discussed 
and  are  con^ared  with  eo^eriment.  Conclusions  are  given  in  Section 
VII. 

II.  BURNIMG  RATE  OF  A  PROPELLAMT  SANDWICH 

In  order  to  see  how  the  burning  rate  of  a  propellant  can  be  mod¬ 
ified  by  placing  it  adjacent  to  a  substance  with  a  different  thermal 
conductivity,  let  us  consider  the  single  system  illustrated  in  Fig.  2.1. 
Two  different  propellants,  each  of  whose  burning  rates  obey  a  pyrolysis 
law,  are  bonded  together.  The  two  propellants  will  be  identified  by  the 
subscripts  2  and  3j  the  subscript  1  identifies  conditions  in  the  gase¬ 
ous  phase.  Let  us  assume  for  sin^jlicity  that  the  specific  heats  of  the 
solids  are  constant  and  equal  (qj  =  c^  =  c  =»  constant),  that  the  densi¬ 
ties  of  the  solids  are  constant  and  equal  (pa  =  Pj  =  p  =  constant),  and 
that  the  thermal  conductivities  of  the  solids  are  constant  (lo  “  constant, 
>,=  constant).  Since  we  are  interested  in  the  effect  of  differing  ther¬ 
mal  conductivities,  we  shall  assume  that  Xg  ^  X,.  The  pyrolysis  law 
governing  the  burning  rate  of  each  propellant  will  be  written  as 
■^a  =  ■'^a  (^1,3)  and  Vg  =  v,  (Tj^g),  where  Vg  and  Vg  are  the  normal  burn¬ 
ing  velocities  (the  linear  regression  rates)  and  Tj  is  the  ten^erature 
at  the  solid-gas  interface.  We  shall  seek  steady- state  solutions  for 
which  the  motion  is  in  the  x  direction  in  Fig.  2.1,  the  temperature  in 
each  solid  depends  only  upon  the  coordinate  nonnal  to  the  surface  of  the 
solid,  and  the  surface  of  each  solid  forms  an  oblique  plane  (as  illus¬ 
trated  in  Fig.  2.1). 

If  V  is  the  velocity  of  the  solid  in  the  x  direction,  then  the 
nonnal  conponent  of  velocity  at  the  surface  is  Vg  =»  v  cos  9g  and 
Vg  =  V  cos  Sg  for  solids  2  and  3«  Hence, 

cos  COS  =  Vg/Vg.  (2.1) 

The  solution  to  the  one-dimensional  heat  conduction  equation  in  a 
moving  medium  has  been  given  in  the  i^endix.  Identifying  the  coordin¬ 
ates  *1  and  f  in  the  i^pendix  with  oxir  present  coordinates  x  and  y,  we 
find  that  £q.  {A.k)  becomes 

Mg  ( X  COB  Gb  y  sin  6(2 ) 

=  Ai  +  Eg  e  (2.2) 
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Figurt  2.1  MODEL  OF  TWO  PROPELLANTS  BONDED  TOGETHER 


Figure  2.2  SANDWICH  PROPELLANT  CONSTRUCTION 
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and 


Ha(x  cos  03  +  y  sin  03) 

T  =  Ag  +  B3  e 
a 


in  phases  2  and  3j  where 

lie  =  p  V  cos  0a  c/Xg  and  p  s  p  v  cos  0  c/Xg. 

3  3 


(2.3) 


Our  objective  is  to  find  constants  and  Bj  and  angles  6^  and 

03  which  satisfy  Eq.  (2.1)  and  the  thermal  conditions  that  must  be 
pliad  at  X  =  -  «>  and  along  the  dividing  plane  y  =  0. 

As  X  the  ten^ierature  in  each  solid  must  approach  the 

initial  teaperature  of  the  propellant,  T*.  Hence, 

4  =  Ag  =  To 

in  Eqs.  (2.2)  and  (2.3).  Furthermore,  the  tea^eratures  must  be  the 
same  in  each  solid  along  the  plane  y  =■  0.  ISiis  can  be  true  at  x  =>  0 
only  if  the  interface  ten^eratures  are  the  same  (Tj^a  “  =  Tj); 

and  it  can  be  true  for  x  <  0  only  if 

^  X  cos  @a  1X3  X  COB  09 

e  =  Ba  e  (2.4) 

for  all  X  <  0.  Equation  (2.4)  will  hold  for  all  x  <  0  only  if 


Efe  =  B3  =  Tj  -  T,  €Uid  COB  4  =  Hg  cos  03. 

In  view  of  the  above  relations  for  |x  ,  Va  and  v  ,  this  last  condi- 

9  9 

tion  lii{)lles  that 


V3^(T,)/v3®(Ti)  =  Xa/Xg, 


(2.5) 


which  is  a  restriction  on  the  physical  properties  of  the  propellants. 
Propellant  pairs  which  do  not  obey  Eq.  (2.5)  may  be  expected  to  have 
more  coii?)licated  ( non-one-dimensional)  temperature  fields.  For  sim¬ 
plicity,  we  shall  consider  only  propellants  for  which  Eq.  (2.5)  is 
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satisfied.  Thus,  Eqs.  (2.2)  and  (2.3)  become 

Pe  cos  ^  ( X  +  y  tan  Gg ) 

Ta  »  To  +  (T,  -  T.)  e 


and 


Me  COB  ^  (x  +  y  tan  6  ) 
T3  =  To  +  (T,  -  T,)  e  * 


(2.6) 


(2.7) 


The  final  thential  boundary  condition  that  the  ten5)erature  field 
must  satisfy  at  the  plane  y  =  0  is  continuity  of  the  normal  heat  flux 
( since  there  are  no  energy  sources  located  at  y  =  O) .  lliis  condition 
can  be  written  as 


(Way)3|y  ^  0  -  (X3T/ay)a|y  ^  0  . 


In  view  of  Eqs,  (2.6)  and  (2.7) »  the  condition  reduces  to 

Xe  tan  Qg  =  X3  tan  63  .  (2.8) 


Equations  (2.1)  and  (2.8)  provide  two  ind^endent  relations  between  the 
unknown  angles  %  and  63,  thereby  determining  each  of  these  angles. 

With  a  little  algebra^  Eqs.  (2.1),  (2.5),  and  (2.8)  can  be  shown  to 
yield 

vAa  “  l/cos  9g  ■  Vl  +  (Xj/^s)  (2.9) 


and 


v/v,  =  l/cos  ©3  =  ^1  +  (>b/^)  . 


(2.10) 


We  have  thus  obtained  unique  angles  0j  and  ©3  such  that  a  siii5)le 
solution  to  the  heat  equation  in  the  solid  phases  exists.  From  Eqs. 
(2.9)  and  (2.10)  it  will  be  noted  that,  if  Xa  >  X^,  then  6g  >  45®  and 
©3  <  45° .  Also,  the  conposite  velocity  v  is  greater  than  both  con^)onent 
normal  burning  velocities  Vg  and  va.  Imbedding  propellant  3  with  a 
high  thermal  conductivity  Xg  in  propellant  2  with  lower  thermal  conduc¬ 
tivity  Xg  thus  increases  the  burning  rate  of  the  original  propellant  2. 
Such  a  sandwich-type  imbedding  is  illustrated  in  Fig.  2.2. 
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III.  MODEL  OF  RATE-ENHANCED  BURNING 

A  reasonable  model  of  rate-enhanced  burning  that  is#  alge¬ 
braically#  somewhat  complicated  has  been  given  in  Ref.  1.  Qualita¬ 
tive  agreement  between  theory  and  experiment  was  obtained.  In  order 
to  clarify  the  essential  physical  characteristics  of  the  process  and 
to  obtain  more  transparent  results#  it  seems  to  be  desirable  to 
analyze  a  model  that  is  much  sinqpler  geometrically  but  retains  the 
physical  phenomena  included  in  Ref.  1.  lliis  model  is  illustrated  in 
Pig.  3.1. 

In  Fig.  3.1#  region  1  is  the  gaseous  phase#  region  2  is  the 
solid  propellant  and  region  3  is  an  inert  material  (bonded  to  the  pro¬ 
pellant)  with  a  high  thermal  conductivity.  A  steady-state#  two- 
dimensional  system  in  which  all  velocities  are  in  the  +x  direction  ia 
inplied.  The  continuity  equation  is  satisfied  automatically  in  this 
model  by  the  statement  that  the  velocity  is  a  constant#  v,  in  the 
solid  regions  (2  and  3)  and  the  velocity  is  Vft|/pi  in  the  gas 
(p  s  density).  IHierefore#  only  the  heat  equations  remain  to  be  con¬ 
sidered. 

Hie  arrows  (  ♦)  in  Fig.  3.1  indicate  the  directions  in  \dilch  heat 
is  supposed  to  be  conducted  in  the  system.  An  adiabatic  (  zero  heat  flow) 
boundary  condition  is  assumed  to  be  applied  tdong  the  plane  y  =  0.  Hie 
inert  solid  3  la  assumed  to  be  thin  enough  and  of  high  enough  thermal 
conductivity  that  the  variation  in  tenperature  across  the  material  (in 
the  y  direction)  is  negligible;  thus#  one-dimensional  heat  conduction 
(in  the  x  direction)#  with  heat  addition  from  the  gas  for  x  >  0  and  heat 
rejection  to  the  propellant  for  x  <  0#  is  assumed  to  occur  in  region  3. 

Hie  gaseous  reaction  zone  will  be  assumed  to  be  so  thin  that  conditions 
in  region  1  can  be  assumed  to  be  practically  unifonn  and  the  tenperature 
T^  can  be  assumed  to  be  approximately  equal  to  the  adiabatic  flame  tem¬ 
perature  T|^.  Actually,  the  highei  velocity  in  the  gas#  caused  by  the  lower 
density  of  the  gas  (see  Fig.  3.2),  leads  to  boundary  layer  development 
(indicated  by  the  dashed  lines  in  region  1  in  Fig.  3.I)  on  the  solid  sur¬ 
faces.  Hie  burning  surface  between  regions  1  and  2  will  be  assumed  to  be  a 
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Figure  3.1  MODEL  OF  RATE  -  ENHANCED  BURNING 


Figure  3.2  VELOCITY  VECTOR  DIAGRAMS  IN  THE  GAS 
AND  IN  THE  SOLID  PROPELLANT. 
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plane.  Therefore,  the  normal  linear  regression  rate  of  the  surface 
is  constant.  Since  we  are  considering  propellants  obeying  a  pyroly¬ 
sis  law,  this  indues  that  the  surface  is  an  isotherm;  l.e.,  the 
surface  ten5)erature  Tj  is  constant.  The  ten5)eratures  in  both  regions 
2  and  3  approach  the  initial  tenperature  To  as  x  At  suffi¬ 

ciently  large  values  of  y,  the  solid  propellant  will  experience  one¬ 
dimensional  heat  conduction  in  a  direction  normal  to  its  surface,  and 
the  ten^erature  profile  will  be  unaffected  by  the  presence  of  material 
3;  [the  only  influence  is  the  si^eriinposed  transverse  velocity  (par¬ 
allel  to  the  surface)  Illustrated  in  Fig.  3>2].  We  shall  assume  that 
this  condition  holds  throughout  region  2a  (divided  from  region  2b  by 
the  dashed  line  normal  to  the  propellant  surface  in  Fig.  3»l)» 
energy  required  to  heat  up  region  2b  is  assumed  to  be  svpplled  by 
heat  conduction  from  material  3« 

The  model  in  Fig.  3*1  will  elLso  be  applicable  to  axially- 
symmetrical  systems  (e.g.,  when  x  is  an  axis  of  symmetry,  region  3 
is  a  cylindrical  wire  and  the  prc^ellant  surface  is  conicea)  provided 
the  thickness  of  the  thermal  boundary  layers  (i.e.,  the  distances  re¬ 
quired  for  T  to  become  practically  equal  to  either  T^  or  To)  in  regions 
1  Euid  2  are  small  conpared  to  the  radius  of  region  3;  numerical  esti¬ 
mates  indicate  that  this  is  often  true  for  wires  greater  than  .001 
inches  In  diameter. 

Consideration  of  the  terperature  field  and  of  the  heat  transfer 
processes  should  deteimine  the  angle  6  in  Fig.  3.1.  IMS  provides  the, 
fractional  Increase  in  the  burning  velocity,  since 

v/vg  “  l/cos  9  (3»l) 

(see  Fig.  3*2;  Vg  s  normal  burning  velocity  of  the  propellant  2). 

IV.  SIMPLE  TIMPERATURE  MATCHING 

In  order  to  obtain  a  quick  and  sinple  (although  inaccurate)  in¬ 
dication  of  how  the  heat  flow  processes  ceui  deteimine  0,  let  us  enploy 
the  following  spproxlmate  Intuitive  criterion:  Along  the  plane  divid¬ 
ing  regions  2  and  3«  the  tenperature  in  medium  3  must  drop  to  T,  in 
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approximately  the  same  distance  that  the  teniperature  in  medium  2  drops 
to  To .  A  gross  violation  of  this  rule  would  lead  to  inordinately 
large  rates  of  heat  transfer  between  materials  2  and  3  in  a  region 
where  the  ten^erature  of  one  of  these  materials  is  practically  uniform. 

An  estimate  of  the  distance  required  for  the  tenq)erature  to 
approach  Tg  may  be  obtained  from  Eq.  (A.3)>  from  which  it  can  be  seen 
that  a  characteristic  decay  length  (in  the  direction  in  which  the  tem¬ 
perature  varies)  for  the  tenperature  profile  in  a  moving  med'.um  is 
l/u.  In  the  propellant,  the  value  of  jx  (in  the  direction  noimal  to 
the  surface)  is 

Ms  =  Ra  Vg  Cg/Xg  =  Pg  V  cos  e  Cg/Xg,  (Jt-.l) 

(c  s  specific  heat  capacity  and  X  s  thermal  conductivity).  Since 
cos  0  =  0  in  the  inert  material  3  (i.e.,  the  teniperature  is  unlfonn  in 
planes  normal  to  the  x  axis). 

Mg  =  Pa  Vg  Cg/Xg  =  Ps  V  Cg/Xg.  (4.2) 

The  geometrical  construction  shown  in  Fig.  4.1  liqplles  that  the  temper¬ 
ature  matching  condition  is 


(l/u  )  cos  0  a  I/mq.  (4.3) 

Equations  (4.1)  through  (4.3)  yield 

l/cos  0  c3/>b)/(p3  Cg/Xg)  .  (4.4) 

The  principeJi  result  exhibited  by  Eq.  (4.4)  is  the  relationship 

l/cos  ®  ~  >  (^"5) 

the  fractional  increase  in  the  burning  rate  is  proportional  to  the 
square  root  of  the  thermal  conductivity  of  material  3.  Essentially  this 
same  result  is  obtained  with  practically  every  reasonable  model  of  the 
process,  and  rough  ejqjerimental  confirmation  of  this  result  has  been 
obtained.  However,  Eq.  (4.4)  does  not  predict  the  observed  d^endence 
of  l/cos  0  iqjon  the  thickness  d/2  of  medium  3-  An  in5>roved  version  of 
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Figur«4.l  GEOMETRY  OF  SIMPLE  TEMPERATURE  MATCHING 
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the  teti5)erature-inatchlng  theory,  given  belov,  does,  however,  yield  a 
thickness  dependence. 

Equation  (A. 3)  really  isn't  applicable  to  medium  3  because, 
as  indicated  in  Fig.  3.1,  heat  is  being  lost  from  medium  3  through 
heat  conduction  to  medium  2.  The  tenperature  of  medium  3  therefore 
drops  to  T,  In  a  distance  that  is  shorter  than  l/viai  with  Hg  given 
by  Eq.  (4.2).  An  estimate  of  the  decay  length,  accounting  for  this 
heat  loss,  may  be  obtained  from  Eq.  (A.7). 

Let  us  assume  that  the  rate  of  heat  loss  is  proportioned,  to 
T  -  Tj .  Then  T^  =  T^  in  Eq.  (A.7).  Since  |i/  >  0,  k'  >  0  and  T  -•  To 
as  X  we  must  have  b  =  0  in  region  3  for  x  <  0.  Hence,  in  region 

[l  +  v^l  +  «  ^  ]  |i'x/2 

3  for  X  <  0,  T  -  To  ~  e  according  to  Eq.  (A.7). 

ISiis  inplies  that  the  characteristic  decay  length  in  medium  3  is 

(2/^/)  [l+yiTT/]  , 

where  and  «/,  defined  in  the  i^pendix,  refer  to  region  3*  Using 
these  definitions,  we  find  that  our  previous  temperature  matching 
condition  becomes 

-1 

[2/(p3  v,  C3/l8)][l  Wi+  (W>^)/(Pa  V,  C3/X3)®]  cos  0 

“  1/(P&  Va  Cb/Xb)  . 

Since  Vj  =  v  =  Vg/cos  6,  where  V3  is  the  normal  linear  regression  rate 
of  the  propellant,  this  formula  can  be  written  as 

cos=  0  =  1  ^/)^)  j  1  Vl  +  XgAfx,  v,  c,)^]cos®  0  j  (4.6) 

where 

'fKX3/(p3  V3  C3)®  =  kVcos®  0 

is  independent  of  0.  Equation  (4.6)  can  easily  be  solved  explicitly  for 
cos®  0;  the  result  yields 
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l/c08  0 


jp2  Og/Xg 

K  Xg 

1  4* 

jpa 

Oa  Cg  ft,  Pg  Vg® 

n-  V2 


Equation  (4.7)  is  the  same  as  Eq.  (4.4)  except  for  the  heat-losa 
correction  term  appearing  inside  the  square  brackets.  OTiis  heat  loss 
is  seen  to  reduce  the  fractional  increase  in  the  burning  rate  (l/cos  0). 
Furthermore,  the  correction  term  has  no  explicit  dependence  upon  1^;  the 
Xg  dependence  given  by  Eq.  (4.5)  remains  valid.  In.  order  to  show  that 
the  correction  term  depends  upon  the  thickness  of  region  3»  ve  must  ob¬ 
tain  an  estimate  of  the  heat  loss  constant  K. 


Since  the  heat  loss  from  material  3  occurs  by  conduction  Into 
material  2,  the  energy  per  unit  sitrface  area  per  second  lost  from  siedlum 
3  is  approximately 

Xe  (T  -  T.)//, 


where  i  Is  a  representative  heat  conduction  length  In  medium  2.  Sie 


energy  lost  per  unit  volume  per  second  is  therefore 


K(T  -  Tj 


h  -  0 

11 

^surface  aree^ 

\  1  i 

M 

^  volume  / 

[>b(T  -  T,)//] 
[X,(T  -  T,)//] 


(2/d)  for 

plane  geometry 
(d)  for  cylin¬ 
drical  geom^ry 


pit 

(4/d) 


Hence, 


K=(4XB)/(/^fjd),  (4.8) 

where  may  be  termed  an  effective  heat  conduction  length  and  d  Is 

either  the  diameter  or  the  thickness  of  region  3>  depending  \Q>on  whether 
region  3  is  a  cylinder  or  a  slab.  As  a  very  rough  estimate,  we  may  set 
-^eff  decay  length  in  region  2,  viz.,  l/Pat  thereby  obtain¬ 

ing 

K  4  Cg  pij  Vg/d  (4.9) 
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[see  Eq.  (4.1)]  and  reducing  Eq.  (4.7)  to 
l/coE  Q  aa  41 -  |1  + 


'Pa 


4 


T  1/2, 


c  Pg  Va  d 


(4.10) 


It  is  clear  from.  Eq.  (4.8)  that  K  depends  upon  the  thickness  d; 
viz.j  K  ~  l/d.  Increasing  d  decreases  K  and  therefore  increases  the 
fractional  enhancement  in  the  burning  rate  l/cos  9.  The  quantity 
l/cos  6  approaches  its  maximum  value,  given  by  Eq.  (4.4),  as  d  -•  a>. 

The  present  formulae  therefore  do  not  predict  the  observed'^  existence 
of  an  optimum  d  and  the  observed  decrease  in  l/cos  0  with  increasing  d 
above  this  optimm  thickness.  The  formulae  developed  in  the  following 
section  do  exhibit  these  refinements.  Equation  (4.10)  predicts  that 
the  dependence  of  l/cos  0  upon  d  enters  through  the  product  Va  d.  This 
result  is  eilso  predicted  by  the  theory  developed  in  the  following  section 
and  is  roughly  confirmed  by  e;q>eriment.^ 

V.  PHENOMENOLOGICAL  THEORY 

In  order  to  develop  a  treatment  that  might  possibly  be  termed  a 
"theory",  let  us  replace  the  teni>eratiu:e  matching  condition  by  an  energy 
balance.  Suppose  that  the  ten^perature  profiles  (along  a  line  normal  to 
the  propellant  surface)  in  region  2  a  of  Fig.  3*1  differ  negligibly  from 
those  in  a  propellant  burning  normally.  Then  no  heat  is  transferred 
across  the  dashed  line  separating  regions  2  a  and  2  b  in  Fig.  3.1.  All 
of  the  energy  required  to  heat  up  the  material  in  region  2  b  must  there¬ 
fore  be  supplied  by  conduction  from  region  3.  Heat  enters  region  3 
only  through  heat  transfer  from  the  gas  in  the  region  x  >  0.  An  over¬ 
all  energy  balance  for  region  3  therefore  states  that  the  energy  entering 
region  3  in  x  >  0  must  equal  the  sum  of  the  energies  required  to  heat  up 
regions  2  b  and  3* 

Let  us  consider  a  slab  geometry  of  unit  height  (normal  to  the 
paper  in  Fig.  3.1) •  The  energy  per  second  required  to  heat  v®  region  3 
then  becomes 
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^  CONHDENTIAL 

^  fii  Vj  d/s  =  Pa  c^  (T^  -  To)  v,  cos  0),(5*1) 


where  the  adiabatic  flame  teii5>erature  is  the  temperature  approached 
in  medium  3  as  x  -*  ®. 

A  close  estimate  of  the  energy  required  to  heat  region  2  b 
may  be  obtained  by  assuming  that  the  temperature  profile  along  the 
dashed  line  bounding  region  2  b  differs  little  from  that  in  region  2  a 
[i.e.,  is  given  by  Eq.  (A.3)]-  If  C  is  the  coordinate  normal  to  the 
propellant  surface  in  Fig.  3*1>  then  (since  T  =  Tj,  the  surface  tempera¬ 
ture,  at  (  a  0  and  T  -*  T,  as  (  -►  -  »)  Eq.  (A.3)  becomes 

Mi  ( 

T  =  T.  +  (T,  -  Tj  e  .  (5.2) 


The  excess  energy  flux  in  region  2bisps  vcg  (T-T, ),  and  the  com¬ 
ponent  of  this  flux  vector  normal  to  the  dashed  line  in  Fig.  3«I  Is 
ftg  V  ca  (T  -  T.)  sin  Q.  The  energy  per  second  required  to  heat  xp  the 
materisil  in  region  2  b  is  therefore 


^  -J" (normal  flux)  d(area)  ”/  Fb 

•  00 

bdry 


V  C8  (T  -  T.)  Bin  0  d(.  (5.3) 


Substituting  Eq.  (5.2)  into  Eq.  (5.3)  and  integrating  yields 


Bb  =  Xb  (Ti  -  T,)  tan  0  ,  (5.4) 

where  use  has  been  made  of  Eq.  (4.1)  for  Mb • 

The  energy  per  second  added  to  material  3  from  the  gas  can  be 
related  to  the  appropriate  heat  transfer  coefficient  hj  the  energy 
per  unit  surface  area  per  second  added  to  medium  3  is  h  (T^  -  T), 
where  T^  is  the  tenperature  of  the  gas  and  T  is  the  local  tenperature 
in  medium  3.  The  tenperature  distribution  in  region  3  for  x  >  0  can 
be  obtained  from  Eq.  (A.J).  Here  T^  =  T^,  and,  since  T  -♦  T^  as 
X  ->  ®,  we  must  have  a  =  0.  If  denotes  the  value  of  T  in  region  3 
at  X  =  0,  then  Eq.  (A.7)  beccmes 
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[l  -  +  »(/]m/x/2 

T  “  -  (T„  -  T^)  e 


(5.5) 


The  quantities  and  «  *  (which  here  refer  to  material  3)  are  defined  in 
the  ^^pendix;  the  value  of  K  sqjpearing  in  the  definition  of  «  ^  is  seen 
geometrically  to  be  related  to  the  heat  transfer  coefficient  throng 
the  fomula 

K=h(2/d)  (5.6) 


for  our  present  slab  geometry.  The  totsQ.  energy  per  second  added  to 
region  3  is 


(5.7) 


where  T  is  given  by  Eq.  (5.5)«  Assvjming  that  h  is  independent  of  x, 
we  find  by  evaluating  the  integral  in  Eq.  (5.6)  that 


El  -  h(T„  -  T  )  (2/p/)  -  1] 


-  1 


2X  cos  9 

h(T  -  T„)  —2 - 

Z  C9  P*  ^3 


^-1 


8h  X,  cos®  9 

1  +  - - - II 


dC  p  Vb 

8  S  “ 


.(5.8) 


The  energy  balance  stated  previously  may  be  written  as 


.  (5.9) 

where  Ej  is  given  by  Eq.  (5«8),  Sg  is  given  by  Eq.  (5-M*  and  E^  is 
given  by  Eq.  (5-l)»  Besides  the  angle  0,  the  only  unknown  parameter 
appearing  in  any  of  these  equations  is  the  tenperature  T„. 

In  order  to  relate  1'  to  9,  we  may  enploy  the  condition  that  the 
energy  to  heat  ipi  region  2  b  is  obtained  by  thennsLl  conduction  from 
region  3*  An  estimate  of  the  amount  of  heat  conducted  cut  of  medium  3 
in  the  region  x  <  0  can  be  obtained  in  teims  of  a  heat  transfer  coeffi¬ 
cient  h/,  defined  in  such  a  way  thart;  the  energy  per  unit  anrea  per  second 
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leaving  medium  3  is  h'  (T  -  T^),  where  T  is  the  ten5)erature  in  medium 
3.  If  we  assume  that  h#  is  Independent  of  x,  then  Eq.  (A.7)  with 
Tj^  =  To  determines  T(x)  in  medium  3  for  x  <  0.  Since  T  =  !„  at 
X  ■  -  oo  and  T  =•  at  X  =  0,  we  find  that  b  =  0  and  a  =  (T^  -  T,) . 

nSierefore  the  total  energy  per  second  leaving  medium  3  is 


h/  (T^  -  To)  e 


[l  +  ,/r+77]  ji/x/2 


dx  =  h/  (T_  -  To)  (2/h/)[i+7FI7] 


-1 


Since  this  heat  flow  rate  must  be  equal  to  Eb>  we  obtain  (en;)loylng  the 
definitions  of  and  k'  and  the  relation  between  K  and  h^) 


2X  cos  6 


i 


1  +  A 1  + 


8h/X-  cos®  o' 


B  -  a  „  a 
dCg  Pa 


-  1 


(5.10) 


An  estimate  of  h^  is  h'  =  le/i  where  £/  is  a  representative  heat 
conduction  length  in  region  2  b.  It  may  be  reasonable  to  assme  that 
I*  is  proportional  to  the  distance  indicated  in  Fig.  If.l;  hence,  let 
us  set/  /  *=  e  sin  0  ''^a  <%)»  where  the  coefficient  e  is  assuned 

to  be  Independent  of  all  experimentally  variable  pctrameters.  Equation 
(5.10)  then  reduces  to 


% 


2Xj  cos  0 
c  sin  0 


f  8  Pb  Cg  ^  cos"  0  ’ 

l  +  Jl+  - 

^  d  Cj®  pj®  Vj  €  Bln  0 


-1 

.  (5.11) 


Equations  (5.l)x  (5.*^)i  (5.8)*  (5.9)  and  (5.11)  serve  to  detezmlne  0. 

\ 

If  we  introduce  Eq.  (5.4)  into  Eq.  (5.11)*  solve  for  (T^  -  T^)* 
and  substitute  the  result  into  Eq.  (5.8)*  we  find  that 


h  (T„  -  To) 


El 


2Xg  cos  0 

=3  Pa  Va 


Vf  8hX  cos®  0 

1  +  — * - 

dC,®  P,®  Vg’ 


/Ti-ToVc, 

fb/Xa/  2 


tan®0  . 


1  + 


/  8rb  ci 

dc,»  p. 


Xj  cos  0 


Vg  e  sin  0 


(5.12) 
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Let  US  introduce  the  dimensionless  variables 

Of  =  (P&  Ca/>B)/(Pa  (5.13) 

[a  measure  of 

y  =  (Ps  ^a'^/(Pa  Ob)>  (5.14) 

^  measure  of 

8  =  d  PS  Va  cs/Xb»  (5.15) 

[a  measure  of  d], 

^  =  (T,  -  Tj/(T^  -  T.)  (5.16) 

and 

H  2  h  d/Xa  .  (5.17) 


When  Eqs.  (5.1)*  (5.*f)  and  (5.12)  are  substituted  into  Bq,  (5.9)>  we 
find,  in  terms  of  these  dimensionless  variables. 


1  -  g  tan»0 

1  +  Jl+  Qp!  ■ 

K  67c  sin  6 

=.  L/rTUZiZi  .1]  . 

L*  6^7  J  L  2  ®. 


(5.18) 


Equation  (5.18)  is  e;q)ected  to  determine  9  in  terms  of  known  parameters. 

In  order  to  use  Eq.  (5.I8)  to  conq)ute  0,  it  is  necessary  to  know 
the  heat  transfer  coefficient  h  so  that  H  may  be  related  to  0,  B,  etc. 
In  the  plane  two-dimensional  system  presently  under  consideration,  beat 
is  transferred  from  the  gas  to  medium  3  8y  convection  through  a  flat- 
plate  boundary  layer.  The  quantity  h  is  then  proportional  to  the  ratio 
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of  Xj  (the  theimal  conductivity  of  the  gas)  to  the  thickness  of  the 
boundary  layer.  The  boundary  layer  thickness  is,  in  turn,  proportional 
to  x/^/^  >  where  Re  =  (v^  -  v)  x/’li  is  the  local  Reynolds  number, 

(Pl  =  density  of  the  gas,  =  velocity  of  the  gas,  ■  coefficient  of 
viscosity  for  the  gas) .  Since  mass  conservation  shows  that  v/vj  ■ 
Pi/pb,  we  find 


Uie  appropriate  average  value  for  x  in  this  expression  is 

(2/n/)  [VI  +  ^  according  to  Eq.  ( 5 « 5) •  Substituting  the 

above  e^qpression  into  Eq.  (5«17)  and  utilizing  Eq.  (^.1^)  yields 

1  -  (Pi/Ph)!^^ 


H  6 


\  1  -  (Pi/fb)! 
ni  ca®  (b  Va  x)  corT~J 


(5.19) 


which  shows  that  the  8  and  9  dependence  of  H  is  determined  by  H  ~  8 
(x  cos  0)"  The  ejqplicit  dependence  of  H  vpon  cos  0,  6,  etc.,  can 

be  unraveled  by  substituting  the  previously  given  expression  for  the 
average  value  of  x,  viz.. 


X  f)i 


2X 

_ 8 

Cs 


cos  0 

p,  ^3 


8h  a  cos° 

8*  7 


(5.20) 


into  Eq.  (5*19)  and  solving  for  H. 

Of  greater  practical  interest  is  the  cyllndrlcally  synmetrical  . 
system.  It  can  be  shown  that,  provided  the  thickness  of  the  thermal 
boundary  layer  in  the  propellant  (region  2)  is  small  coupared  with  the 
diameter  of  region  3j  the  only  modifications  in  our  previous  develcp- 
ment  are  that  an  addltioned.  factor  of  l/2  appears  in  the  expression  for 
E^  [on  the  right-hand  side  of  Eq.  (5.1)]  and  that  an  additional  factor 
of  2  ‘ppears  in  the  numerator  on  the  right-hand  side  of  Eq.  (5.6).  This 
last  modification  leads  to  an  additioneJ.  factor  of  2  in  the  nume..ator 
Inside  the  radical  sign  in  Eqs.  (5.8)  and  (5.II),  thus  causing  the  co¬ 
efficient  8  to  be  rqxlaced  by  16  inside  both  rsuUcal  signs  in  Eq.  (5.I8). 
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Qlie  first  modification  causes  the  factor  1/2  to  be  replaced  by  ijh  in 
the  last  teirn  in  Eq,.  (5  *18). 

The  ejqjression  for  h  (and  therefore  for  H)  in  the  cylindrically 
symmetrical  case  is  the  same  as  that  given  above  for  the  plane  case, 
provided  the  thickness  of  the  gaseous  boundary  layer  on  medium  3  is 
small  coii5)ared  with  the  diameter  of  medium  3  thro'oghout  the  region  in 
which  appreciable  heat  transfer  from  medium  1  to  medium  3  occurs. 
Nmerical  estimates  indicate  that  this  condition  will  be  satisfied  for 
d  ~  a  few  thousandths  of  an  inch  provided  x  [as  estimated  from 
Eq.  (5.20)]  is  less  than  roughly  1  cm.  Since  a  direct  estimate  of  x 
from  Eq.  (5*20)  yields  x  ~  a  few  cm.,  the  validity  of  the  above  condi¬ 
tion  is  questionable,  and  it  may  be  better  to  en^ploy  eaipirical  formulae 

2 

for  h  for  small  wires  placed  parallel  to  the  flow  direction.  Jacob 
quotes  resxilts  of  this  type  obtained  by  Muhler  in  turbulent  flow  with 
the  Reynolds  number  based  upon  wire  diameter,  Re^,  in  our  present  range 
of  interest  (to  <  Re^  <  toOO);  these  results  have  also  been  corrobor¬ 
ated  (in  a  close  cq)proximation)  by  some  very  recent  work.  Qhe  resvilts 

quoted  by  Jacob  may  be  correlated  roughly  by  the  formula  Nusselt 

JaT 

nvonber  =  hd/Xj  sa  .3  Re^  ,  a  reasonable  and  convenient  eq>proxdjnation 
to  which  is 

.5 

h  =  ,25  (Xa/d)[pi(vi  -  v)  d/nj  .  (5.21) 

Since  pj  <  <  r,,  we  have  v  <  <  Vj,  and  Eq.  (5 .21)  can  be  written  as 

H  =  p  6’Vcos’^  e, 

where 

®  ‘  -"5  fe)  (^) 

and  use  has  been  made  of  Eqs.  (5*15)  and  (5.17)  and  of  the  identity 
Pj  v^  =  Pjj  V  =  pfe  Vg/cos  0.  Although  the  accuracy  inherent  in 
Eq.  (5.22)  is  not  great,  this  relation  does  provide  a  usefxil  order- 
of -magnitude  estimate,  and  it  appears  to  be  as  accurate  an  expression 

-  19  - 


(5.22) 

(5.23) 


CONFIDENTIAL 


CONFIDENTIAL 


1 
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as  one  can  obtain  for  a  flow  system  in  which  t,he  /jooinetry  is  as  com¬ 
plicated  as  that  presently  under  consideration. 

By  substituting  Eq.  (5*22)  into  Eq.  (5«l8)  and  also  Introducing 
the  other  modifications  (indicated  above)  that  arise  in  Eq.  (5*l8) 
for  cylindrical  geometry,  we  obtstin 


-  [V-  -  ^  • 


(5.2h) 


Solution  of  Eq.  (5*24)  (by  trial  and  error)  for  the  value  of  9  corres¬ 
ponding  to  prescribed  values  of  a,  P,  7,  &,  e  and  t  determines  the 
fractional  enhancement  in  the  bviming  rate.  In  subsequent  sections  we 
shall  consider  only  cylindrically  symmetrics^L  systems,  enploylng  Eq. 
(5*24)  instead  of  Eq.  (5.18). 

VI.  RESULTS  AND  COMPARISON  WITH  EXPERIMENT 

One  very  general  consequence  of  Eq.  (5*24)  is  that  the  dependence 
of  9  upon  both  d  and  Vj  enters  only  through  the  parameter  6j  i.e.,  a, 

P,  7,  e  and  r  are  all  independent  of  both  d  and  Vj  [see  Eqs.  (5*13) > 
(5.14),  (5*16)  and  (5"23)3.  This  result  actually  does  not  even  dq>end 
upon  the  heat  transfer  assun5>tions  associated  with  the  cylindrical 
geometry  since  Eqs.  (5-18)  through  (5*20)  can  be  shown  to  yield  the  same 
result  for  plane  gecsnetry.  In  view  of  the  definition  of  6  [Sq.  (5»15)], 
we  may  conclude  that  0  is  a  function  of  the  product  dva  when  all  other 
parameters  are  held  fixed.  The  normal  linear  regression  rate  Vg  can  be 
varied  practically  Independently  of  all  other  psurameters  appearing  in 
our  formulae  by  varying  the  chamber  pressure.  E)q)erimental  data  on  the 
effect  of  varying  both  d  and  Vg  independently  is  therefore  available. 

A  plot  of  l/cos  0  vs.  dvg  for  various  values  of  Vg ,  obtained  from  the 
ejqjerimental  data,  is  shown  in  Fig.  6.1,  A  snail  systematic  influence 
of  Vg  is  noticeable,  thus  indicating  that  the  foundations  of  our  theory 
are  not  entirely  sppllcable  to  these  experiments.  In  extreme  cases 
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Figurt  6.1  -  The  ratio  of  the  enhonced  burning  rote  to  the  normal  burning  rote  os  o 
function  of  the  product  of  the  normol  burning  rate  ond  the  diometer, 
computed  from  the  experimental  data  of  reference  1. 
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(near  the  peak  value  of  l/cos  6),  the  discrepancy  is  seen  to  amount  to 
nearly  a  factor  of  2  in  l/cos  Q  for  given  dVg .  Tnis  same  objection  may 
be  raised  against  the  more  con^licated  theory  of  Ref.  1.  There  are  a 
number  of  possible  causes  of  the  discrepancy;  for  example,  we  did  not 
consider  the  fact  that  if  the  number  of  wires  and  the  chamber  size 
remain  fixed  as  the  wire  diameter  is  increased,  then  the  over-all 
adiabatic  flame  temperature  decreases  as  d  increases.  Our  primary  rea¬ 
son  for  emphasizing  this  discrepancy  is  to  illustrate  that  the  present 
theory  cannot  reasonably  be  expected  to  yield  enhanced  burning  rates 
that  will  agree  with  ordinary  ejperiments  more  closely  than  a  factor  of  2. 

Various  general  characteristics  of  Eq.  (5*24)  are  of  interest. 

For  exanple,  it  can  be  seen  that,  as  6  -►0,  the  factor  in  the  curly 
brackets  in  Eq.  (5*24)  must  vanish,  and  Eq.  (5«24)  therefore  yields 

v2/3 


1 

cos  6 


1  => 


1  /«  7  &\ 

2  Wr-  ) 


(6.1) 


for  very  small  values  of  6.  It  can  also  be  shown  from  Eq.  (5*24)  that, 
as  8  -•  •,  cos  0  again  approaches  unity;  explicitly. 


cos 


-  1  a 


.-3/2 


(6.2) 


for  very  large  values  of  6.  To  obtain  Eq.  (6.2)  involves  retaining  a 
three-teim  expansion  of  the  radical  on  the  right-hand  side  of  Eq.  (5-24). 

Since  — i— »  -»  1  as  8  -•  0  and  as  6  we  conclude  that  l/cos  0  attains 
cos  0  J  -t 

at  least  one  maximum  value  for  positive  6.  The  quantity  l/cos  0  is 
bounded  for  all  8^0  because  the  quantity  in  the  curly  brackets  in 
Eq.  (3.24)  must  remain  non-negative;  setting  the  curly  bracket  equal  to 
zero  yields 

cos®  0-1  (1  “  ®) 


(6.3) 


(providing  a  cubic  equation  for  the  minimum  value  of  cos*  0),  thus 
implying  that 
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l/cos  6  ^  1  (6A) 

for  all  6^0  [since  the  right-hand  side  of  Eq.  (6.3)  is  non-negative]. 
The  maximum  value  of  l/cos  0  must  therefore  be  finite.  Equation  (6.4) 
shows  that  the  bound  on  the  fractional  increase  in  the  burning  rate  is 
proportional  to  [i«e.,  proportional  to  »/  1^/^  ;  see  Eq.  (5-13)] 
for  the  physically  in5)ort€int  range  of  values  of  the  parameters  (viz., 
for  o/r€  >  10^)j  this  result  has  already  been  suggested  in  Section  IV. 
[see  Eq.  (4.5)]  and,  as  has  previously  been  indicated,  is  in  rough  agree¬ 
ment  with  experiment. 

jBie  value  of  l/cos  0  as  a  fionction  of  &,  obtained  from  Eq.  (5*24) 

3  -1 

for  the  representative  case  a  =  10  ,  P  =  10  ,  7  =  1,  c  =  1,  r  =  1,  is 

plotted  in  Fig.  6.2.  The  values  of  5  in  the  experiments  reported  in 
Ref.  1  lie  between  1  and  50.  Conparison  of  Figs.  6.1  and  6.2  indicates 
qualitative  agreement  between  theory  and  experiment.  IHie  order  of  mag¬ 
nitude  of  the  maximum  enhancement  in  the  biuming  rate  is  the  same  on 
the  theoretical  and  experimental  curves,  llie  largest  difference  be¬ 
tween  theory  and  experiment  appears  to  occur  at  small  values  of  6  and 
can  be  reduced  considerably  by  suitably  adjusting  the  values  of  c 
and  r  • 

VII.  CONCLUSIONS 

The  model  given  in  Section  III  and  the  analysis  given  in  Section 
V  appear  to  expladn  the  underlying  physiceJ.  mechanism  of  rate-enhanced 
burning.  The  principal  result  of  this  study  is  Eq.  (5.24),  which 
determines  the  bixming  rate  as  a  function  of  the  physiceO.  parameters  of 
the  problem.  Uhe  primary  qualitative  conclusion  is  that  the  ratio  of 
the  enhanced  burning  rate  to  the  noimaJ.  bximing  rate  is  prcportionail  to 
the  square  root  of  the  ratio  of  the  thermal  conductivity  of  the  wire  to 
that  of  the  prcpellant  (i.e.,  to  ./o’  ).  For  ordinary  materials,  this 
square  root  does  not  exceed  about  50,  and,  when  the  correct  constant  of 
proportioneiLity  is  enployed,  it  is  found  that  the  burning  rate  can  be 
Increased  at  most  by  a  factor  of  10.  It  may  also  be  noted  [see^  for 
exanple,  Eq.  (6.4)]that  a  small  value  of  r  (i.e.,  a  small  ratio  of  the 
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Figure  6.2  -  The  ratio  of  the  enhanced  burning  rate  to  the  normal  burning  rote  oe  a 
function  of  the  dintensionless  diometer  8,  comfMted  from  eq.  (5.24),  for 
a«IO*,  y»l,  *■!,  T»l. 
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propellant  surface  temperature  to  the  adiabatic  flame  ten^jerature) 
favors  a  larger  relative  increase  in  the  burning  rate. 

There  are  many  related  basic  research  problems  that  to 

be  interesting.  For  exan^le,  an  analysis  of  steady-state,  reacting, 
flat-plate  boundary  layer  flow,  with  a  solid  phase  located  upstream 
(at  the  forward  edge  of  the  plate)  and  with  the  plate  moving  slowly  in 
the  flow  direction,  woxild  appear  to  be  tractable  (although  perh£5)s 
difficult)  and  would  represent  the  flow  in  the  vertex  where  the  pro¬ 
pellant  meets  the  wire.  However,  no  analyses  of  this  type  are  essential 
for  predicting  the  degree  of  enhancement  of  the  burning  rate;  the  con¬ 
trolling  physical  mechanism  of  rate-enhanced  burning  Is  now  well- 
understood. 
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APPENDIX 

SOLUTIONS  OF  THE  STEADY-STATE  ONE-DIMENSIONAL  HEAT  CONDUCTION  BODATION 

IN  A  MOVING  SOLID 


Let  us  consider  a  material  for  which  the  density  p,  the  heat 
capacity  c  and  the  thermal  conductivity  X  are  all  constant.  If  this 
material  is  moving  with  velocity  v,  then,  in  the  absence  of  heat  sources 
or  sinks,  the  steady-state  differential  equation  governing  the  distribu¬ 
tion  of  the  tenperature  T  in  the  material  becomes 

cpv*VT-XV®T=0.  (A.l) 

For  a  non-rotating  solid  body,  v  is  constant  everywhere;  i.e.. 


V  -  V  e^  , 

where  v  is  the  (constant)  magnitude  of  the  velocity  and  e^  is  a  unit 
vector  in  the  V  direction.  Let  us  investigate  systems  in  which  the 
ten^)erature  is  a  function  only  of  a  single  cairtesian  coordinate  {  . 

If  e^  denotes  a  unit  vector  in  the  (  direction,  then  Eq.  (A.l)  reduces 
to 


(c  p  v/x)  (e^  •  e^)  dT/d  f-  d®T/d{®  »  0  . 


(A.2) 


Ihe  dot  product  egf^earlng  in  Eq.  (A.2)  is  simply 


where  Q  is  the  angle  between  the  (  and  directions.  It  may  be  noted 
that  the  coiic>onent  of  v  in  the  (  direction  is 


e^  “  v  cos  0  =  Vj 
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We  shall  assume  that  cos  d  /  0  and  shall  define  a  characteristic  length 


ii  =  p  V  cos  9  c/x. 

The  general  solution  to  Eq.  (A.2)  then  becomes 

T  =  A  +  B  e^^^,  (A.3) 

where  A  and  B  are  arbitrary  constants.  In  the  plane  containing  e^  and 
e^  we  may  introduce  the  cartesian  coordinate  system  illustrated  in 
Fig.  A.l;  [cj.  is  normal  to  e^  and  t  =  *i  cos  0  +  f  sin  s].  In  terms  of 
the  coordinates  i|  and  t,  Eg..  (A.3)  la 

T  -  At  B  e‘^(’'  cos  0  +  f  sin  0)  ^ 

Ohls  result  is  used  in  Section  II  and  elsewhere. 

Next,  let  us  remove  the  restriction  that  no  heat  sources  or 
sinks  can  be  present  but  in5)ose  the  additional  condition  that  e,j  and  e^ 
are  parallel  (i.e.,  0  =  0  and  the  {  and  axes  coincide).  Thus,  we 
shall  consider  one-dimensional  steady-state  heat  conduction  with  heat 
addition  or  extraction  and  with  motion  of  the  solid  in  the  direction  of 
the  tengjerature  gradient.  We  shall  assume  that  the  energy  per  unit 
volume  per  second  added  to  or  extracted  from  the  solid  is  given  by 
K  (Tj^  -  T),  where  K  (energy/vol.  sec.  is  a  positive  constant  and 
T„  is  a  constant  reference  ten?)erature.  Energy  is  added  to  the  solid 

X\ 

if  T  <  Tg,  while  energy  is  extracted  from  the  solid  if  T  >  T^.  The 
one-dimensional,  steady-state  heat  equation  now  becomes 

c  p  V  dT/dx  -  Xd®  T/dx®  =  K  (Tj^  -  T),  (A.5) 

^ere  x  is  the  coordinate  in  the  v  direction.  In  teims  of  the  reciprocal 
of  the  characteristic  length. 
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Figure  A. I  THE  ADOPTED  CARTESIAN  COORDINATE  SYSTEM 
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Ji/  =  p  V  c/x 

and  the  dimensionless  heat  transfer  parameter^ 

=  (4k/x)/h/2  , 

Eq.  (a. 5)  can  be  vritten  in  the  form 

p/d  (T  -  Tp)/dx  -  d®  (T  -  Tj^)/dx®  +  (p/®K/A)  (T  -  Tj^)  -  0. 
The  general  solution  to  Eq.  (A. 6)  is 

T  -  T  +  a  ef^  ^  xOp^S  +  ^  "  /T+~7f  ]p'::/i' 

R 

where  a  and  b  are  art>ltrary  constants.  Equation  (A.7)  is  used  in 
Sections  IV  and  V. 
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